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We present results from a comprehensive analytical and numerical study 
of nonequilibrium dynamics in the 2-dimensional complex Ginzburg-Landau 
(CGL) equation. In particular, we use spiral defects to characterize the do- 
main growth law and the evolution morphology. An asymptotic analysis of 
the single-spiral correlation function shows a sequence of singularities - analo- 
gous to those seen for time-dependent Ginzburg-Landau (TDGL) models with 
0(n) symmetry, where n is even. 
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Much recent interest has focused on pattern formation in the complex Ginzburg-Landau 
(CGL) equation: 

^^^^ = + {l + ^a)V'^P{r,t) - (1 + ^/3)|V^(r, t)|2^(f, t), (1) 

where il){r^t) is a complex order-parameter field which depends on space (r*) and time (t). 
In Eq. (|lD, a and j3 are real parameters. The CGL equation arises in a range of diverse 
contexts, as reviewed by Cross and Hohenberg This universality arises from the fact 
that the CGL equation provides a generic description of oscillations in a spatially-extended 
system near a Hopf bifurcation p|. 

The CGL equation exhibits rich dynamical behavior with variation of the parameters 
a and /3, and the "phase diagram" has been investigated (mostly numerically) in various 
studies 1^. In a large range of parameter space, the emergence and interaction of spiral 
defect structures characterizes the morphology. In this letter, we study the nonequilibrium 
dynamics of the CGL equation resulting from a small-amplitude random initial condition. 
In general, this nonequilibrium evolution is referred to as "phase ordering dynamics" or 
"domain growth" , and constitutes a well-studied example of far-from-equilibrium statistical 
physics 1^,^ . Our analytical understanding of phase ordering systems has depended critically 
upon modeling the dynamics of defects in these systems (e.g., interfaces, vortices, monopoles, 
etc.) §-[11. In this letter, we use spiral defect structures to characterize the evolution 
morphology in the CGL equation. Many novel features emerge in our study, which should 
be of great relevance for both experiments and subsequent numerical simulations. 

For simplicity, we will focus on the CGL equation with a = and dimensionality d = 2. 
However, the results presented here are also relevant for the cases with a ^ and d > 2, as 
the underlying paradigm does not change, i.e., spirals continue to determine the morphology 
in large regions of parameter space. Following the work of Hagan [^j, Aranson et al. ||10|| , 
and Chate and Manneville 0], let us briefly discuss the phase diagram of the d = 2 CGL 
equation with a = 0. The limiting case (3 = corresponds to the dynamical XY model, 
which is well understood. The appropriate (point) defects are vortices, and domain growth 
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is driven by the attraction and annihilation of vortex-antivortex pairs. The relevant growth 
law for the characteristic length scale is L{t) ~ (t/lnt)^/^ W^M i the analytic form 
of the time-dependent correlation function (which characterizes the evolving morphology) 
has been obtained by Bray and Puri, and (independently) Toyoki ||^. Without loss of 
generality, we focus on the case with /? > 0. For < /? < /?i (/?i ~ 1.397 p), spirals 
(which are asymptotically plane-waves) are linearly stable to fluctuations. For j3i < (3 < (32 
{P2 — 1-82 0,11), spirals are linearly unstable to fluctuations, but the growing fluctuations 
are advected away, i.e., the spiral structure is globally stable. Finally, for P2 < P, the spirals 
are globally unstable and cannot exist for extended times [T^. Our results correspond to 
the parameter regime with (3 < (32- 

Figure 1 shows the typical evolution from a small-amplitude random initial condition 
for the case with (3 = 0.75. Our numerical simulations were performed by implementing an 
isotropic Euler-discretization of Eq. (|l|) on A^^-lattices (A^ = 256 for Figure 1), with periodic 
boundary conditions in both directions. The discretization mesh sizes were At = 0.01 and 
Ax = 1.0. In Figure 1, we plot constant-phase regions and the relevant color-coding is 
provided in the figure caption. The evolving morphology is characterized by spirals and 
antispirals, and there is a typical length scale L, e.g., inter-spiral spacing or the square root 
of inverse defect density, which is the definition we will use subsequently. 

Figure 2(a) plots ln[L(t)] vs. Int for 5 representative values of (3. The length-scale 
data was obtained from 5 independent runs on A^^-lattices with = 1024. After an initial 
transient period, the length scale L{t) should saturate to an equilibrium value (Lg) because 
of an effective spiral-antispiral repulsive potential [|l|. This should be contrasted with the 
(3 = case, where vortices continue to anneal (at zero temperature) as t ^ 00. As a matter 
of fact, the data for (3 = 0.25, 0.5 in Figure 2(a) does not exhibit this morphological freezing 
on the time-scales of our simulation, though signs of the crossover are evident for (3 = 0.5. 

To understand this crossover behavior, we recall the analytical solution for an m-armed 
spiral due to Hagan M: 
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■^(r, t) = p{r) exp [—iut + im6 — i4>{r)] , (2) 

where r = (r, 6); and = — q"^), where g is a constant which is determined by (3 0. The 
hmiting forms of the functions p(r) and 0(r) are 

p(r) ^Jl — q^, (t^\^) Qi as r oo, 

p(r) — i> ar™, — r, as r ^ 0, (3) 

where a is a constant which is determined by finiteness conditions. We will focus on the case 
with m = ±1, as only 1-armed spirals are stable in the evolution 0. Furthermore, we are 
only interested in distances r ^ ^, where ^ is the defect core. Thus, we consider the spiral 
form in Eq. (|^) with p(r) = y/1 — q'^ and 0(r) = qr (appropriate for r oo). 

We expect that spirals behave similarly to vortices for L < L^, where qL^ ~ 0(1). Thus, 
the early evolution should be analogous to that for the XY model, both in terms of the 
domain growth law and correlation function. In Figure 2(a), the solid line has a slope of 1/2 
and the initial growth (at least for (3 < 0.75) appears to be consistent with the behavior for 
the XY model, i.e., L(t) ~ (t/ Int)^/^ for d = 2. We also expect the saturation length Lg to 
scale with Lc- Figure 2(b) plots Lg vs. q^^ for a range of /3-values, and demonstrates that 
our numerical data is consistent with Lg ~ q~^. We can also obtain the scaling law for the 
crossover time and the corresponding numerical results (not shown here) are in agreement 
with it. 

Next, we consider the correlation function for the evolution morphology shown in Figure 
1. It is obviously relevant to first consider the correlation function for a single spiral of length 
L, as the snapshots in Figure 1 can be thought of as consisting of disjoint spirals of size 
L. (Of course, this ignores modulations of the order parameter at spiral-spiral boundaries 
but we will discuss those later.) We have approximated the 1-armed single-spiral solution as 
ip{f, t) ~ y/1 — q'^ exp [—iut + i{6 — qr)]. The correlation function is obtained by considering 
the correlation between points fi and f2 (= ri + fi2) and integrating over fi as follows: 

C{ru) = ^J dmeWn,t)^{r2,ty}h{L-r2) 
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Re / dfi exp [i{9i -62- qri + g|fi + fi2|)] h{L - |fi + fi2|) 



(4) 



V 



where V is the spiral volume; and we have introduced the step function h{x) = 1 (0) if 
a; > (a; < 0). The step function ensures that we do not include points which lie outside 
the defect of size L. 

It is convenient to introduce variables 9i — 9^ = 9; x = ri/L; r = ri2/L, to obtain 



where we have used V = nL^ in d = 2. Thus, the scaling form of the single-spiral correlation 
function is C(ri2)/C(0) = g{ri2/L,q'^L'^). In general, there is no scaling with the spiral size 
because of the additional factor qL. We recover scaling only in the limit q = {P = 0), 
which corresponds to the case of a vortex. Essentially, spirals of different sizes are not 
morphologically equivalent because there is more rotation in the phase as one goes out 
further from the core. 

Figure 3 plots C(ri2)/C(0) vs. ru/L for the case with (3 = 0.75 [q ~ 0.203). These 
results are obtained by a direct numerical integration of Eq. (^. We consider 4 different 
values of L. The functional form in Figure 3 exhibits near-monotonic behavior for small 
values of L (i.e., in the vortex or XY limit); and pronounced oscillatory behavior for larger 
values of L, as is expected from the integral expression. Notice that ri2/L < 2 - larger 
values of correspond to the point r2 lying outside the defect. 

The asymptotic behavior of the correlation function in the limit r = ri2/ L —>■ (though 
f'u/C, ^ 1) is of considerable importance as it determines the tail of the momentum-space 
structure factor [^. In particular, we are interested in the singular part of the correlation 
function as r ^ 0. In this limit, we can discard the step function in Eq. (^ as it only 
provides corrections at the edge of the defect. The asymptotic analysis of the integral in Eq. 
(^) involves considerable algebra, which we will report in detail elsewhere. Here, we confine 
ourselves to quoting the final result for the singular part of C(ri2): 




exp —iqLlx — {x'^ + r'^ + 2xrcos9y^'^\ h[l ~ {x'^ + r'^ + 2xr 008 9)^^"^] 



(5) 
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2p=o^o (2p)!(2m)!r(i _p) (^ + ^ + 1)12 

(2m + 1) {2p + l)r2(™+P+i) In r. (6) 

Eq. (H) is one of the central results of this paper and we would like to briefly discuss its 
implications. The leading-order singularity is the same as that for the XY model {jS = q = 0), 
Csing('"i2) = |r^lnr [T^, as expected. However, there is also a sequence of sub-dominant 
singularities proportional to (gL)^r'^lnr, (gL)'^r^lnr, etc., and these become increasingly 
important as the length scale L increases. These sub-dominant terms in Csmgi'f'u) are 
reminiscent of the leading-order singularities in models with 0{n) symmetry, where n is 
even 0,0. Of course, in the context of 0{n) models, these singularities only arise for n < d 



as there are no topological defects unless this condition is satisfied. In the present context, 
all these terms are already present for d = 2. The implication for the structure-factor tail 
is a sequence of power-law decays with S{k) ~ (qLy^"^^^^ L'^ /(kL)'^'^'^"^, where m = 1,2, 
etc. Thus, though the true asymptotic behavior in (i = 2 is still the generalized Porod tail, 
S{k) ~ L'^{kL)~'^, it may be difficult to disentangle this from other power-law decays. 

Finally, Figure 4 compares our numerical data for the correlation function with the 
functional form of the single-spiral correlation function. Recall that the correlation function 
does not scale with the characteristic length because of the spiral nature of the defects. In 
Figures 4(a)-(c), we have plotted numerical data for C(ri2, t)/C(0, t) vs. ri2 at t = 500, and 
(3 = 0.75, 1.0, 1.25. For the comparison with Eq. (^, the length scale L is taken to be an 
adjustable parameter. In each case, the best-fit value of L matches the length scale obtained 
from the inverse defect density (see Figure 2(a)) within 10 percent. As is seen from Figure 
4, the single-spiral correlation function is in good agreement with the numerical data for the 
multi-spiral morphology upto (approximately) the first minimum. As a matter of fact, the 
agreement is excellent (perhaps fortuitously) for (3 = 1.25, shown in Figure 4(c). 

In the context of phase ordering dynamics, the Gaussian auxiliary field (GAF) ansatz 
P]- [S] has proven particularly useful for the characterization of multi-defect morphologies. 
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We have critically examined the utility of the GAF ansatz in the present context |T3[ and 
find that it is only reasonable at early times - where, in any case, the ordering process is 
analogous to that for the XY model. We are presently studying methods of improving the 
GAF ansatz for the CGL equation and will discuss this elsewhere. 

More generally, the utility of the GAF ansatz arises from the summation over phases 
from many defects, which results in a near-Gaussian distribution for the auxiliary field. 
However, in the present context, the shocks between spirals effectively isolate one spiral 
region from the influence of other regions. As a matter of fact, the waves from other spirals 
decay exponentially through the shock and the phase of a point is always dominated by the 
nearest spiral. Therefore, we expect that the correlation function will be dominated by the 
single-spiral result - in accordance with our numerical results. 

To summarize: we have undertaken a detailed analytical and numerical study of nonequi- 
librium dynamics in the CGL equation. For early times (L < Lc ~ q~^), the domain growth 
process is analogous to that for the XY model, which is well understood. At later times, 
distinct effects due to spirals are seen and the evolving system freezes (in a statistical sense) 
into a multi-spiral morphology. We have undertaken an asymptotic analysis of the correla- 
tion function C(ri2) for a single spiral. It exhibits a sequence of singularities as ri2/L 0. 
Furthermore, this correlation function is in good agreement with the numerical data for 
multi-spiral morphologies, over an extended range of distances. 
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FIGURE CAPTIONS 



Figure 1: Evolution of tlie CGL equation from a small- amplitude random initial condi- 
tion. The evolution pictures were obtained from an Euler-discretized version of Eq. (|l]) 
with a = 0, (3 = 0.75, implemented on an A^^-lattice (A^ = 256). The discretization mesh 
sizes were At = 0.01, Ax = 1.0; and periodic boundary conditions were imposed in both 
directions. The snapshots show regions of constant phase 6^ = ta.n^^{lmilj/Iieip), measured 
in radians, with the following color coding: 9^ G [1.85,2.15] (black); 6^ G [3.85,4.15] (red); 
9^ G [5.85,6.15] (green). The snapshots are labeled by the appropriate evolution times. 

Figure 2: (a) Plot of \n[L{t)] vs. Int for a = and p = 0.25,0.5,0.75,1.0,1.25 - de- 
noted by the specified symbols. The characteristic length scale, L{t), is obtained from the 
square root of the inverse defect density - measured directly from snapshots as shown in 
Figure 1. The numerical data shown here was obtained as an average over 5 independent 
runs for A^^-lattices (with = 1024). The solid line has a slope of 1/2. 
(b) Plot of saturation length Lg vs. q^^ for a range of /5-values. The corresponding values 
of q (as a function of (3) are obtained from Hagan's solution, cf. Figure 5 of Ref. ||^. The 
solid line denotes the best linear fit to the numerical data. 

Figure 3: Correlation function for the 1-armed spiral solution when (3 = 0.75 (g ~ 0.203). 
We plot C(ri2)/C(0) vs. ri2/L for different spiral sizes, L = 15, 25, 50, 100 - denoted by the 
specified line-types. These results are obtained from a direct numerical integration of Eq. 
(I)- 

Figure 4: Numerical data for the correlation function C{ri2,t)/C{0,t) vs. ri2 at t = 500 
for the cases a = and (a) P = 0.75; (b) /5 = 1.0; (c) [3 = 1.25. The numerical data 
was obtained as an average over 5 independent runs for A^^-lattices (with N = 1024). The 
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solid line refers to the numerical integration of Eq. (|^) with L as an adjustable parameter. 
Subsequently, the ri2-axis is scaled so that the point C(ri2, t)/C(0, t) = 1/2 is matched for 
the numerical data and the analytical expression. 
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This figure "figl.gif" is available in "gif" format from: 



http://arXiv.org/ps/cond-mat/0103527vl 



